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Abstract 

In order to explain and model the inner ring in photoluminescence (PL) patterns of indirect excitons in 
GaAs/AlGaAs quantum wells (QWs), we develop a microscopic approach formulated in terms of coupled 
nonlinear equations for the diffusion, thermalization and optical decay of the particles. The origin of the 
inner ring is unambiguously identified: it is due to cooling of indirect excitons in their propagation from 
the excitation spot. We infer that in our high-quality structures the in-plane diffusion coefficient D^ 2d) - 
10 - 30cm 2 /s and the amplitude of the disorder potential Uo/2 0.45 meV. 

PACS numbers: 73.63.Hs, 78.67.De, 05.30.Jp 
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More than two decades ago, long-distance diffusion and drift transport of charge-neutral exci- 
tons with a long lifetime was optically visualized in bulk Si and Cu 2 [LlL l2U - hi these earlier works, 
strain gradient potential traps were used to induce the drift motion. Due to the low particle con- 
centrations, the transport was described in terms of a classical picture, with no quantum-statistical 
corrections, and the diffusion coefficient was attributed to exciton-phonon scattering. In-plane 
propagation of long-lived indirect excitons in coupled QWs over large distances has also been 
reported iQ, 0, Q, U Bl : in this case one has a unique possibility to optically map the quasi-two- 
dimensional motion of composite bosons. Furthermore, the density n 2t j of indirect excitons can be 
large enough to ensure nonclassical population of the ground-energy state, Ne=o = e To/T - 1 > 1, 
where T = {n% 2 n 2 d)l{2k B M x ) and T are the degeneracy temperature and exciton temperature, re- 
spectively, and M x is the exciton in-plane translational mass (for a recent review of cold exciton 
gases in coupled QWs see Ref. Jst). 

One of the most striking features of photoluminescence associated with indirect excitons is the 
appearance of two PL rings Jo]. While the second, external ring has already been explained in 
terms of in-plane spatially separated electrons and holes 1*1 [ill, the origin of the first, inner ring, 
which arises purely due to the transport of indirect excitons, remained unclear. The underlying 
physical picture we propose to explain this ring is that in the optically-pumped area the exciton 
temperature T is much larger than the lattice temperature T b . As a result, the optical decay of 
excitons is suppressed, but while they diffuse out they cool down and eventually become optically- 
active, giving rise to a local increase of the PL signal. Thus the first ring is a generic feature of cw 
PL from excitons, direct or indirect, that are nonresonantly excited in high-quality QWs. 

In this Letter we develop a microscopic theory for the long-range transport, thermalization and 
optical decay of QW excitons, model the inner PL ring, show an effective screening of QW dis- 
order for n 2 d > 10 10 cm -2 , and analyze the quantum-statistical corrections. By numerically fitting 
the PL spectra we clarify the main scattering channels which contribute to the diffusion of indirect 
excitons and evaluate the diffusion coefficient and amplitude of the QW disorder potential. In the 
proposed model, the exciton temperature T is affected by heating due to the optical excitation, 
heating due to the LA-phonon assisted conversion of the mean-field energy into the internal en- 
ergy, and recombination heating or cooling due to the optical decay of low-energy excitons. The 
mean-field energy of indirect excitons also gives rise to a potential energy gradient and, therefore, 
to the in-plane drift motion. 

In the experiments, the PL pattern is imaged by a nitrogen-cooled CCD camera with the spectral 
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selection by an interference filter adjusted to the indirect exciton energy. Fine adjustment of the 
filtering energy is done by rotating two interference filters relative to the optical axis. As a result, 
we are able to remove the low-energy bulk emission that otherwise dominates the spectrum near 
the excitation spot 0, and observe a two-dimensional spatial image of the inner ring (see Figs, la, 
lc, and Id). Note that the inner PL ring can be missed if the bulk emission is not removed from the 
PL signal. In Figs, le and If we plot the measured exciton PL in the energy-coordinate plane. The 
exciton energy £p L decreases with increasing distance from the excitation spot, as detailed below. 
This results in an arrow-shaped profile of the exciton PL images in the E-x coordinates. The 
external ring is also seen at high excitations, both in x-y (Fig. la) and E-x (Fig. le) coordinates. 
The excitation is done by a HeNe laser at 633 nm (the laser spot is a Gaussian with FWHM 
= 2<t - 6 /jm, the excitation power P ex = 1 - 400 //W). The coupled QW structure with two 8 nm 
GaAs QWs separated by a 4 nm Alo.33Gao.67 As barrier is grown by molecular beam epitaxy. 

Our approach to the transport, relaxation and PL dynamics of indirect excitons is formulated in 
terms of three coupled nonlinear equations: a quantum diffusion equation for n 2 d, a thermalization 
equation for T, and an equation for the optical lifetime T opt = l/T opt . The quantum-statistical 
corrections, which enhance the nonlinear effects, are included in the description. 

The nonlinear quantum diffusion equation is given by 

dn 7 



'2d 



= v[z^ 2d) Vn 2d + jufVdV (u n 2d + U QW ) ] - Y opt n 2d + A , (1) 



dt 

where r opt , D^ A \ fi®\ an d A are the radiative decay rate, diffusion coefficient, mobility, and 
generation rate of QW excitons, respectively, Uq W = t/ r and( r ||) is a random potential due to the QW 
thickness and alloy fluctuations, and the operator V is defined in terms of the in-plane coordinate 
vector ry. The generalized Einstein relationship Jill . //^ = D^ d) [(e To,T - l)/(k B T Q )], which 
yields the classical limit f/ 2d) = D^ d) /(k B T) for T » T , strongly enhances the nonlinearity of the 
diffusion Eq. © for T < T . The positive mean-field energy u n 2d on the right-hand side (r.h.s.) 
of Eq.© is due to the well-defined dipole-dipole repulsive interaction of indirect excitons 112 . 
[l^lllJ]- Here uq = 4n(e 2 /sb)d z , £b is the background dielectric constant, and d z is the separation 
between electron and hole layers. The mean-field energy gives rise to the in-plane drift motion 
with the velocity = -jr x ) a Vn 2 d- As a result, an effective screening of the disorder potential 
t^rand by dipole-dipole interacting indirect excitons builds up with increasing n 2d llli I15I : the 
excitons tend to accumulate near the minima of C/ r and( r ||) [local increase of Mo^2d( r ii)] an d avoid the 



maxima of £/ ran d(r||) [local decrease of Mo"2d(r||)]. As we show below, in our high-quality structures 



U (Q) = 2<|[/ rand (ru)|) « 0.9 meV and the mean-free energy utfifl - 1.6 meV for nt® = 10 10 crrT 2 , 
so that at low exciton temperatures T ~ 1 K the QW disorder is strongly screened and practically 
removed for n 2A Z n 2 }. 

The temperature dynamics of excitons is given by 



d _ I dT \ S pump + S op t + S d 
dt ~ \dt) m + 2k B Th - k B T I 2 ' 
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where 7i = (l-e- To/r ) Jz[z/(e z + e~ ro/r -l)], 7 2 = e" r ° /r Jz[(ze z )/(e z +e- r ° /r -l) 2 ]. Equation 
© describes the LA-phonon assisted thermalization of indirect excitons: r sc = (;r 2 ft 4 p)/(7>jp.M^v s ) 
is the characteristic scattering time, E = 2M x v 2 , v s is the longitudinal sound velocity, p is the 
crystal density, and D dp is the deformation potential of exciton - LA-phonon interaction. The form- 
factor F z (x) = [sin(x)/x][e' x /(l -^ 2 /tt 2 )] refers to a QW confinement potential and characterizes 
a spectral band of bulk LA-phonons which interact with indirect excitons, a = (dQwM x v s )/h, 
and Jq W is the QW thickness. For 1.5 K < T < 6.5 K and n 2 d = 10 10 cm -2 , relevant to our 
experiments, Eq. © yields a thermalization time 0.15 ns > r th > 0.03 ns. The strong increase of 
T th with decreasing T is due to classical slowing down of the relaxation kinetics, which occurs at 
k B T ~ E (excitons with energies E < E /4 cannot emit LA-phonons) Hi. 

Heating of indirect excitons by the laser pulse is given by the term 5 pump = {E- mc -k B TI 2 )K TQ > 
on the r.h.s. ofEq. ©. Here, A To = [(nh 2 )/(2kBM x )]A(t,r\\). The generation of indirect excitons 
is a secondary process, mainly due to quantum tunnelling of photoexcited direct excitons to the 
energetically more favorable states of indirect excitons. The excess energy E inc of a created (in- 
coming) indirect exciton is large: it exceeds the energy splitting between the direct and indirect 
excitons, which is about 20meV. For the highest generation rates used in the experiments, the 
exciton temperature T m3X =*■ 6.4 K at the laser spot centre is much larger than 7b- 1 .5 K. 

Recombination heating or cooling of QW excitons is given by the term 5 opt = \k B TI 2 Y 0Vt - 
EjTop^To on the r.h.s. of Eq. ©, where 

i / e \ r 1 i + z 2 

r oP t = HHH p peiIt i dz (4) 

2r R \ k B T ) Jo Be~ z E y ^ T - 1 
and r^, = [E y /(2r R k B T )] f Q l [(l - z 4 )/(Be- z2 ^ /kBT - 1)] dz are the optical decay rates for the 
concentration and energy density of indirect excitons, respectively, r R is the intrinsic radiative 
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FIG. 1: (a), (c)-(f) Experimental and (b) calculated patterns of the PL signal from indirect excitons. In (a), 
the PL intensity in the area within the green circle is reduced by a constant factor for better visualization. A 
bright spot in the middle of the inner ring in (a) is due to residual bulk emission, (e) and (f) Image of the PL 
signal in the E-x coordinates. The external PL ring is also seen for high excitations, both in (a) x-y and (e) 
E-x coordinates. For (a) and (e) the excitation power is ,P ex = 250 /uW and for (b)-(d) and (f) P ex = 75 /iW, 
respectively, and T\, = 1.5 K. 

lifetime of indirect excitons, E y = ^/(2M X ), B = (e E ? /kBT )/(l - e~ To/T ), and p = (E x ^/s^)/c 
(E x is the total energy of a ground- state indirect exciton). In contrast with the evaporative cooling 
schemes used in atomic optics to remove high-energy atoms from magnetic traps 

Q], the 

optical evaporation of QW excitons is an inherent process, which deals with the lowest-energy 
particles, < E < E y , from the radiative zone. Both signs of 5 opt can be realized: 5 opt > 
(S opt < 0), i.e., recombination heating (cooling) of indirect excitons for k B T > E y and &b7o 
(k B T «; E y and &b7o)- In our case recombination heating occurs: far away from the excitation 
spot T exceeds 7b by about 3 mK, i.e., due to relatively high 7\,, the effect is small. Note that an 
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FIG. 2: The PL intensity /pl, (a) measured and (b) calculated with Eqs. against radius ry for four 

optical excitation powers P ex . The Gaussian profile of the optical excitation, with cr = 3 /mi, is shown by 
the thin solid lines. The cryostat temperature 7\, = 1.5 K. In numerical evaluations the following parameters 
are used: Dd p - 9.6 eV, tr = 13 ns, E- mc = 20meV, M x = 0.215 mo, c?qw = 8nm, and d z = 11.5nm. The 
best fitting parameters are Uq - 0.9 me V, D^f^ - 60cm 2 /s, and C x _ x - 15cm 2 /s. The upper inset: the 
energy position of the PL line, £pL = £pl(/||)- The lower inset: the measured (square points) and calculated 
(solid line) inner ring radius versus P ex . The triangular points refer to the external PL ring. The measured 
(circle points) and calculated (solid line) HWHM spatial extension of the PL signal against r\\. 

effective recombination cooling of particles takes place for T b < 0. 1 K: the effect is strong and 
reduces r th by > 20 %. 

Heating of excitons by in-plane drift is given by 5 d = -7oWo(v to t- Vn 2 d) on the r.h.s. of Eq. ©, 
where v tot = Vdiff + Vdrift and v d iff = -(D x 2d) ln 2 d)Vn 2 d. This heating mechanism is due to conversion 
of the mean-field energy into the internal one. The effect is particularly important for k^T » Eq 
(E /k B x 0.4 K), when the momentum relaxation time is much less than r th . 

In our model, the diffusion coefficient D x 2d) = (Di 2d !D (2d) )/(D? d ! + D m ) has two con- 

' x v x-x x -imp // v x-x x-imp' 

tributions: diffusion due to scattering by imperfections (QW impurities and bulk LA-phonons), 
with D (2d) = Df I (7b), and self-diffusion due to exciton-exciton scattering. The latter 

x imp x imp 
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channel is important for n 2d > 10 10 cm" 2 and is approximated by = C X - X (T/T ) 

Jill . For r\\ far away from the excitation spot, the asymptotic solution of Eqs. 0-© yields 
/ PL oc exp[-(T opt /D^_)) l/2 r\\\. In contrast, the experimental data show a much more steep de- 
cay of the PL signal and its spatial pinning at a critical radius r™ = r"(P cx ) (e.g., 40//m 
for P ex = 402 /iW, see Fig. 2a). We attribute such a behaviour to the n 2 d-dependent screening of 
long-range-correlated QW disorder by dipole-dipole interacting indirect excitons. The narrowing 
effect is illustrated in Fig. 3a for a particular realization of the disorder potential t/ r and(ll)- hi or- 
der to include the long-range-correlated disorder we use a thermionic model [11] which operates 
with n 2d - and T-dependent diffusion coefficient Z) x 2d) = D X A) exp[- U^ 0) /(k B T + u n 2 d)]. Note that 
for r|| - rf, where the narrowing effect is weak due to small n 2 d, a more adequate description 
of the diffusion and thermalization kinetics deals with phonon-assisted hopping between the QW 
localized states 

In order to model the experimental results within the developed microscopic picture, we solve 
Eqs. ©-© numerically for a stationary, cylindrically-symmetric optical excitation profile, so that 
the generation rate A(r||, f) = A(ry) oc P ex exp(-r 2 /cr 2 ). The best fit for the experimental data 
plotted in Fig. 2a yields U Q = 0.9 me V, D^f^ = 60cm 2 /s, and C x _ x = 15cm 2 /s. The calculated 
spatial profile of the PL signal, 7 PL = 7pl(1|), is shown in Fig. 2b for various pump powers P ex . 
While the density profile n 2 d = «2d( r i|) always has a bell-like shape (see Fig. 3b and the inset of 
Fig. 2a), with increasing P ex the inner PL ring develops in the / PL -profile. This is in a complete 
agreement with the observations. The inner ring has a nearly classical origin, and arises due to 
heating of indirect excitons by the optical excitation [5 pump term in Eq. ©]: with increasing r\\ the 
exciton temperature T rapidly decreases towards T b (see Fig. 3b); as a result, the optical lifetime 
T op t = 1 /r op t decreases too (see the inset of Fig. 3b), giving rise to a local increase of /pl0"||) at 
r\\ = r^. Thus the inner ring is a spatial counterpart of the PL-jump observed in the time-resolved 
experiments 112 ill . Our numerical simulations also reproduce the observed increase of rjj g and 
the spatial extension of the PL area (HWHM of the signal) with increasing P ex (see the inset of 
Fig. 2b). 

The finding of the fitting parameters, which refer to the total diffusion coefficient D X A \ is 
complex, i.e., we fit all the curves plotted in Fig. 2a [I PL = 7pl(^|) and E PL = E PL (r\{) for various 
Pex] by using the same values of Uq, D x _), and C x _ x . The blue shift £ PL > of the PL energy 
£pL is due to the mean-field interaction energy of indirect excitons, 6 PL = E PL - E x = u n 2d (see 
Figs, le-lf and the inset of Fig. 2a). Thus we use the measured £ PL to estimate n™*(r\\ = 0), and 
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FIG. 3: (a) The narrowing effect: screening of the long-range-correlated disorder potential C/ r and(r||) by 
dipole-dipole interacting indirect excitons. The QW effective potential C/ e ff(ll) ~ ^rand(ll) + u o n 2d(r\\) 
calculated with Eqs. 0-© for harmonic U ran d(r\\) with Uq/2 = 0.45 meV (see the inset), (b) T = T(r\\) 
(solid line) and «2d - «2d(/||) (dashed line). Inset: T opt = T opt (r||). (c) Diffusion velocity Vdiff = VdiffCll) ( s °lid 
line) and drift velocity Vdrift = Vdnft(li) without (dashed line) and with (dotted line) quantum-statistical 
corrections. Inset: Ne=o = Ne=o(t}\)- Plots (a)-(c) refer to the maximum excitation, «™ ax = 2.5 x 10 10 cm -2 
(see the solid lines in Fig. 2). (d) Diffusion coefficient D^ d) = D^ d \r\\) for n™* = 0.16 x 10 10 cm~ 2 (solid 
line), 0.44 x 10 10 cm -2 (dashed line), and 1.27 x 10 10 cm -2 (dotted line). These values of «™ ax correspond 
to the PL signal shown in Fig. 2. 

therefore the generation rate A(ry = 0), necessary for numerical modelling with Eqs. The 
amplitude of the disorder potential Uq determines the steepness of hiSj\\ > rjj 8 ) and dependence 
r \\ = r y r ^ ex -' °^ tne ^ pinning radius. In turn, D^f) an( * ^ x ~ x determine the ring contrast and 
r \\ = r y g ^ ex ^ dependence. The total diffusion coefficient D^ d) = Dx 2d) (ry) is plotted in Fig. 3d. 

In Fig. 3c we show that for the small size excitation spot (cr - 3 /urn) used in the experiments, 
the drift velocity Vdrift, due to the gradient of the mean-field interaction energy Wo«2d(/||X is much 
larger than the diffusion velocity Vdiff. The total velocity has a maximum value v™ t ax (ry rjj 6 ) - 
1.5 x 10 5 cm/s for n™ x (r\\ = 0) - 2.5 x 10 10 cm"" 2 (see Fig. 3c). Note that in our case the mean- 
field energy gradient uo\Vn 2 d{r\\ - r?)\ ^1.6 eVVcm exceeds the maximum strain-induced gradient 
\VU\ - 0.4eV/cm used in the experiments yj,|2|]. 
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Finally, we emphasize that in our experiments, which deal with the cryostat temperature T b = 
1.5 K, nonclassical occupation numbers of modest values, iV^(r|| r|j g ) =t 1, build up at the 
position of the inner ring (see the inset of Fig. 3c). As illustrated in Fig. 3c, in this case the quantum 
statistical corrections, e.g., to Vdiff and to the Einstein relationship, are about 35% and therefore 
cannot be neglected. Nonclassical statistics occurs at the position of the inner ring, where the 
exciton gas is already cold but still dense. Furthermore, for 7\, ~ 0.1 K (not yet realized in an 
optical imaging experiment) numerical modelling with Eqs. ©-© gives well-developed Bose- 
Einstein statistics with A^Cry r|j g ) » 1. 

We appreciate valuable discussions with L. V. Keldysh, L. S. Levitov, L. Mouchliadis, and B. D. 
Simons. Support of this work by EU RTN Project HPRN-2002-00298 is gratefully acknowledged. 
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